Theory of frequency-dependent spin current noise through correlated quantum dots 
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We analyze the equilibrium and non-equilibrium frequency- dependent spin current noise and spin 
conductance through a quantum dot in the local moment regime. Spin current correlations are 
shown to behave markedly differently from charge correlations: Equilibrium spin cross-correlations 
are suppressed at frequencies below the Kondo scale, and are characterized by a universal function 
that we determine numerically for T — temperature. For asymmetrical quantum dots dynamical 
spin accumulation resonance is found at the Kondo energy, u> ~ Tic- At higher temperatures 
surprising low- frequency anomalies related to overall spin conservation appear. 



PACS numbers: 72.25.-b, 73.63.Kv, 72.15.Qm, 72.70. 

Introduction. Coherent detection and manipula- 
tion of spin currents in nanostructures has recently at- 
tracted wide attention due to possible promising appli- 
cations in future storage technologies and quantum com- 
puting [l|, 0] ■ Many proposals have been made to build 
spin batteries to inject spin-polarized current, and then 
filter, manipulate, and detect it Q. Often one makes use 
of ferromagnetic electrodes in these circuits while in 
other cases the application of an external magnetic field 
5 1 or the presence of a ferromagnetic resonance process 
6J enables one to filter and detect spin currents. Quan- 
tum dots play a special and important role in this regard: 
In these devices, the strong electron-electron interaction 
enables one to manipulate the spin of a single electron Q , 
and such quantum-dot devices provide a possible route 
to quantum computing 

However, to use spin circuits efficiently, it would be of 
crucial importance to characterize the noise in them. In 
addition, the structure of the noise provides valuable in- 
formation on interactions and correlations. In fact, a lot 
of attention has been devoted to noise analysis in corre- 
lated mcsoscopic circuits for this reason (jJ-QjJ. Due to 
progress in experimental technology, it is now possible to 
measure ac conductance properties as well as frequency- 
dependent noise in these circuits down to very low tem- 
peratures, and even in the Kondo regime [12l 4l4l - Fur- 
thermore, with efficient spin filtering methods [3[ mea- 
suring spin-resolved current noise in such circuits is also 
within reach. Surprisingly, while a lot is known about 
the properties of ordinary noise in the correlated regime, 
much less is known about the structure of spin current 
noise. So far, only spin correlations in the sequential tun- 
neling [TBI , and perturbative regimes 17 1 have been 
analyzed, and these works focused almost exclusively on 
shot noise. 

Here we carry out a detailed analysis of the full fre- 
quency spectrum of the spin-dependent current noise in 
the Kondo regime. We show that equilibrium spin cur- 
rent correlations are characterized by two universal func- 



tions, which we determine numerically for T = tem- 
perature usin g th e method of numerical renormalization 
group (NRG) [18( 1 . and compute analytically for large fre- 
quencies. At finite temperatures, we analyze spin corre- 
lations using a perturbative approach. We find in all 
regimes that correlations between electrons of the same 
and opposite spins behave markedly differently. In the 
perturbative regime these remarkable differences emerge 
at frequencies below the Korringa relaxation rate: while a 
dip appears in the frequency-dependent noise of opposite 
spin directions, a large peak develops for the parallel spin 
components. These surprising features are all intimately 
related to spin conservation. 

Model. Focusing on the Kondo regime, we shall as- 
sume that there is a single spin 5 = 1/2 electron on the 
quantum dot, which couples to the electrons on the leads 
through the Kondo interaction [lj| , 
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Here a stands for the three Pauli matrices, the fields 
ip ra = J D D c rr ?{s)de destroy electrons of spin a in leads 
r E {L, R}, and their dynamics are governed by the non- 
interacting Hamiltonian, Hq = ^2 ra J e c\ a (e)c ra {e)de 
[25I ]. The coupling j in Eq. ([T]) is the usual dimension- 
less coupling, which incorporates already the density of 
states in the leads, and is related to the Kondo tempera- 
ture as Tk ~ D e _1 / J , with D the cut-off energy appear- 
ing in f/Vo-- The dimensionless hybridization parameters 
are given by vj, = cos(</>/2) and vr = sin((/>/2), with <fi 
parametrizing the asymmetry of the dot: <f> — tt/2 cor- 
responds to a symmetrical quantum dot with maximum 
transmittance. 

Equilibrium noise. In view of the special structure 
of Eq. (J), it is natural to introduce the 'even' and 'odd' 
linear combinations, = cos(</>/2) ipT,+sm((j)/2) i/jr, and 
^ = sin(^>/2) ipz—cos((f)/2) ipR. Although only ^ couples 
to the spin in H- m t, changing the chemical potential in 
one of the leads couples the fields 'F and <F, and both 
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contribute to the spin noise. 

To compute the noise, we first define the spin compo- 
nent a of the current in lead r through the equation of 
motion, J ra = e N ra = e i[H[ nt , N ra ]. The correspond- 
ing current is found to have two distinct (even and odd) 
parts, J ra = I ra + I r<r , with 



I r a = ejy r i(F^ a -*tF a ) , 
la =ej%i(F^ a -& a F a ) , 



(2) 



and the prefactors defined as ^l/r = [1 ± cos(</>)]/4 and 
1l/r = ±sin(</>)/4. The operator F a = (Sct4') (J denotes 
the so-called composite fermion operator [20I ] , and repre- 
sents the universal (Kondo) part of the dot-electron. 

The operator identity, + I r ± = 0, and the sim- 
ple even-odd decomposition of J ra imply that, in equi- 
librium and in the absence of external magnetic field, 
the sixteen components of the symmetrized noise S°°, = 
\{{Jra{t), Jrv(0)}), depend on just two universal func- 
tions, s and s. Maybe the most interesting left-right noise 
component, S1 R , can be expressed, e.g., as 

e 2 

S lr(u) = -— T K sin 2 (0) (S aa 's(u}) + aa' s(uj)) , 

where e 2 /2ir = e 2 /h denotes the universal conductance 
unit, and the dimensionlcss functions s and s depend 
exclusively on the ratios lo /Tk and T /Tk- The function ,s 
is related to the 'even' current component, and it governs 
the correlations between spin up and spin down carriers, 
however, its contribution cancels in the charge noise and 
charge conductance, which are exclusively determined by 
the 'odd' component of the current, incorporated in s. 

In equilibrium, the fluctuation-dissipation theo- 
rem relates S™, (w) to the real part of the spin- 
conductance through the dot, Re G™ (u>) — 
— — coth(w/2T) (u>) , which can therefore also be 
expressed in terms of two dimensionlcss universal con- 
ductance functions, g(u,T) and g(uj,T). The left-right 
conductance, e.g., reads 

Re Gli (lu) = |- sin 2 (0) (<W$(«, T) + aa' g(u, T)) . 

Using tedious but straightforward manipulations, we 
can express g and g in terms of the spectral functions 
qf{w, T) and Qx a x , (w, T) of the composite fermion and 
of the "current" operator I CT = i{F^„ - F a ), 
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with /(w) denoting the Fermi function |21[. Since F a 
and X a are local operators, we can compute g and g (and 
thus s and s) by using the powerful method of numerical 

rcnormalization group (nrg) dum. 
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Ik- 
fig. 1: (color online) Zero-temperature universal functions s 
and 5 computed by NRG. Inset: universal spin conductance 
functions g and g. 



T = 0, equilibrium results. The T = temperature 
universal functions s{ui/Tk), and s(u/Tk) and the con- 
ductance functions giw/Tx) and g(uj/Tx) are displayed 
in Fig. [T] The high-frequency behavior of s and s can be 
captured by doing perturbation theory in j and summing- 
up the leading logarithmic corrections to give 



s{lo/T k ) m -§ s{u/T K ) 



3-n 2 H 1 
16 T K ln 2 (\u>\/T K ) 



for bj 3> Tk- Though they look similar at high frequen- 
cies, s and s behave markedly differently in the Fermi 
liquid regime, lo <C Tk, where s = ol\uj\/Tk + ■ ■ ■ , while 
s = a (\uj\/Tk) 3 + ■ ■ ■ , with a and a universal constants 
of the order of unity. The lu 3 scaling of s is related to 
spin conservation: In the absence of external spin relax- 
ation mechanism, the total number of spin up electrons 
can fluctuate between two values, Nf and TVf + 1. Since 
the spin up electrons couple to the spin down electrons 
only at a single point (the quantum dot), no steady spin 
current can be generated for the spin down electrons by 
injecting spin up electrons in one of the leads. Thus the 
spin conductance G^, (ui) must vanish at ui = 0, and by 



analyticity, G 



n 



In equilibrium, however, the spin 



current noise is simply related to the spin conductance by 
the fluctuation-dissipation theorem, implying a |u;| 3 scal- 
ing of s at T = 0. This argument carries over to finite 
temperatures too, where it leads to an asymptotic behav- 
ior, s ~ ru 2 in the absence of external spin relaxation. 
We should emphasize that, in our calculations, spin re- 
laxation is due to the interaction part of the Hamiltonian 
which, however, conserves the total spin, and leads to the 
vanishing of spin noise component at lj = 0. Introduc- 
ing some source of an external spin relaxation, however, 
leads to a violation of spin conservation, and amounts in 
a finite = 0)^0 0. 

The fundamental difference between ft an d t4- cor ~ 
relations shows up even more strikingly in the spin- 
conductance (see Fig. QJ: While G^ R (u>) is dominated 
by g(co) and behaves qualitatively the same way as the 
conductance through the dot, G]^ R (u>) ~ g(ui) exhibits a 
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FIG. 2: Diagrams contributing to the noise. The reduced 
density matrix of the spin evolves along the upper and lower 
Keldysh contours. Triangles denote the bare current vertices, 
dots indicate the exchange interaction. Arrows correspond 
to conduction electron propagators. To leading order, the 
current-current correlation function is given by the connected 
diagrams (top), S CO nn(u}), and the "disconnected" diagrams 
(second line), S'disc(w), with II the propagator describing the 
evolution of the reduced density matrix of the spin (third 
line). The dressed current vertex W ra is given by diagrams 
similar to those of the self-energy, £ (last line), with one of 
the dots replaced by a triangle. 



resonance at a frequency u> « 0.5 Tk [26J. This can be 
understood in a simple and intuitive way: The spin con- 
ductance tJ- is generated by flips of the localized spin. 
For u> > Tk, the coupling to the conduction electrons 
gets stronger with decreasing uj, and increases the con- 
ductance. At very small energy scales, uj >C Tk, however, 
the impurity spin is quenched, and with the above mech- 
anism being absent, the t4- conductance must vanish. 

T / 0, perturbative regime. Computation of the 
finite temperature noise requires care: Usual finite tem- 
perature NRG broadening procedures lead to an unphys- 
ical finite linear coefficient for s(w), conflicting with our 
exact finite T result, s(uj) ~ uj 2 . Therefore, for T / 0, 
other methods must be used. For T ^> Tk, we carried 
out a systematic expansion in j for the time-dependence 
of the reduced density matrix of the spin and the spin 
current noise using the formalism of Rcfs. 



241. De- 



tails of this involved calculation shall be published else 
where 
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here we just outline the main results. 
Naively, to calculate the noise in leading order, one 
would just compute the first (connected) noise diagram 
of Fig. [3J S conn (uj). This diagram accounts for short 
time current correlations mediated by electron-hole ex- 
citations in the leads, and dominates indeed the noise 
at high and intermediate frequencies, uj > T. At small 
frequencies, however, a resummation of the perturbation 
series is necessary, because there the "disconnected" con- 
tribution, 5disc(^), turns out to be of the same order in 



j as S conn (uj), and becomes also important: This con- 
tribution accounts for correlations between subsequent 
incoherent tunneling processes, generated by the impu- 
rity spin itself. These correlations are due to the mere 
fact that a spin flip process where a conduction electron's 
spin is flipped from up do down, f— h|,, must be followed 
by a process 4-— >t- To account for them, one needs to 
solve a Dyson equation for the propagator II of the re- 
duced density matrix of the spin, as sketched in Fig. [2] 
In this approach, spin relaxation is characterized by the 
relaxation rate, r(w), appearing in the self-energy £ of 
the propagator IT [231 ], 



T' 



T 



(4) 



where Re L(x,y) — ir y + n[x sh(a;) — y sh(y)]/[ch(x) — 
ch(y)], and Im L(x, y) — ^ / dx'He L(x' , y)/{x — x'). In 
the uj — > limit, T(uj) can be identified as the Korringa 
relaxation rate, Ek = r(0)/2, of the impurity spin, which 
for a simple voltage-biased quantum dot reads 



E 



K 



irj 2 T 



eV , eV 
—= coth — — 
2T 2T 



In the voltage-biased case, we can express the left-right 
component of the spin current noise in a compact form 



S^(o;) = --Re 



oJ_ r 2 M -r 2 (uj) 

16 -iu + T(uS)/2 



(5) 
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with R(uj) = j 2 cos((f>) L(ui/T, 0). The symmetrized noise 
is shown in Fig. [3] At high frequencies, ui 3> Ek, the 
noise is dominated by the result of simple-minded per- 
turbation theory, corresponding to the second line of 
Eq. ([5]). This part of the correlation function describes 
short-time correlations within a single tunneling process, 
generated by the dynamics of electron-hole excitations in 
the leads. However, at time scales t ~ 1/lu > 1/Ek, con- 
secutive incoherent tunneling processes start to correlate 
by the constraint mentioned before. These correlations 
are captured by the first term in Eq. ([5]), coming from 
the "disconnected" part of the noise (see Fig. [2j). As 
a consequence, for uj < Ek, a large dip appears in the 



noise component Sjj^, while a bump emerges in 
zero-bias, V = 0, we find that S^ r (oj = 0, V = 0) = 0, in 
agreement with the fluctuation-dissipation theorem and 
the observation that the linear spin conductance between 
spin up and spin down electrons must vanish. This has 
a simple physical explanation: A spin-t electron injected 
from the left can give rise to a spin-J, outgoing electron on 
the right with a certain probability (see Fig. [3]). However, 
before such a process occurs again, another spin flip pro- 
cess must take place, where the dot spin is flipped back. 
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FIG. 3: (color online) Top: Sketch of consecutive spin flip 
processes. Bottom: Equilibrium and non-equilibrium noise 
spectra (lj) in the perturbative regime, max{T, ui, Ek} > 
Tk, as computed from a diagrammatic approach. In the non- 
equilibrium case a simple current bias was assumed, V£ = V 
and Vr = 0, while j — 0.07 and <f> — n/2 in both cases. 



In equilibrium, this second process (on the average) re- 
moves exactly the same amount of 4- spin from the right 
lead as injected in the first process. Therefore, no equi- 
librium 'hi. dc spin conductance is possible. 

Remarkably, the above correlations only show up in 
the spin current noise, and cancel out in the charge 
current noise, Slr = J2a a' ^lr ■ Being the result of 
rather classical correlations between subsequent incoher- 
ent processes, these low frequency features can also be 
captured by a much simpler rate equation approach (sec 
Refs. [l7, 21 1), which, however, is unable to account for 
the high frequency part of the noise at ui > T . 

Although the above results are perturbative in j, 
they carry over to the whole regime max{T, u, Ek} > 
Tk with the small modification that j must be re- 
placed by the renormalized coupling j —> j(T,ui,eV) « 
l/]n(max{T,uj,E K }/T K ). 

Fermi liquid regime. In the Fermi liquid regime, 
lo, T, eV >C Tk, one can compute the spin current corre- 
lations by describing the dot in terms of scattering states 
that interact at the impurity site. This is a rather cum- 
bersome approach for finite frequencies. However, ob- 
serving that correlations between spin up and down elec- 
trons are generated only through the residual electron- 
electron interaction, simple phase space arguments im- 
mediately give that the T = shot noise is just given 

'k, while for 



ture and finite frequency noise and the precise determi- 
nation of these universal constants is very complicated, 
and shall be considered in a future publication. 

Conclusions. Analyzing the full frequency depen- 
dence of the spin current noise through a quantum dot 
in the Kondo regime we found that ti correlations are 
strongly suppressed at frequencies below the Kondo tem- 
perature and below the Korringa relaxation rate as com- 
pared to fp correlations due to overall spin conserva- 
tion. In the t4- conductance a resonance is predicted 
at uj ~ Tk- Observing these striking features is within 
reach with present-day noise measurement techniques. 
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by Sf R {V) = (e 2 /27r) 7 sin\4>){eVf/T 1 
equilibrium we recover the numerically observed result, 
Sl l R (uj) = (e 2 /2n) a sin 2 (<f))\uj\ 3 /T K , with 7 and a two 
universal numbers. The discussion of the finite tempera- 
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